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We present empirical features of parton energy loss in nucleus-nucleus collisions at RHIC through 
studies of the spectra and nuclear modification factors (Raa) for charged hadrons, neutral pions (tt ) 
and non-photonic electrons. The flat distribution of Raa at high transverse momentum (pr) for a 
given collision centrality is consistent with a scenario where parton energy loss Apr is proportional 
to pr- The centrality dependence of the parton energy loss indicates the absence of path length 
O ' dependence in the magnitude of energy loss. The lack of strong path length dependence suggests 

a dynamical picture where the dense partonic medium undergoes rapid expansion and the density 
f"*^ . of the medium falls rapidly in the first a few Fermi interval, which may be much shorter than the 

full path length. Implications of the empirical constraints on the parton energy loss will also be 
discussed. 

PACS numbers: 25.75.Ld 
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I Nucleus-nucleus collisions at RHIC produce a hot and dense medium. Particle production from these collisions 

. exhibited various distinct physical phenomena depending on transverse momentum (px) scales. Below 2 GeV/c, 
^ ' hydrodynamic calculations can describe the expansion of the strongly interacting matter formed in the collisions 

Between 2 and 5 (or 6) GeV/c particle production exhibits a Constituent Quark Number (CQN) scaling reflecting 
coalescence or recombination hadron formation mechanism with underlying constituent quark degrees of freedom [2j • 
^| ', Particle production in the higher px region is consistent with features of parton fragmentations, and may be used 
■ to study parton energy loss in the hot and dense medium. Experimental measurements at RHIC have shown a 
£SJ \ strong suppression of high px particles with respect to the expectation of binary nucleon-nucleon collision scaling for 
CN i charged hadrons [|[ , for neutral pions 0, Q and for non-photonic electrons @ . Early parton energy loss calculations 
have focused on gluon radiative energy loss mechanism which has been used to explain the high px suppression of 
light hadrons Recent measurements of the strong suppression of non-photonic electrons from heavy quark 

decays contradict the expectation of reduced radiative energy loss for heavy quarks where a dead-cone effect has been 
calculated [9|. It has been proposed that collisional energy loss is important for heavy quarks and cannot be neglected 
[lol. [III. fl2Ul3j . However, a consistent dynamical description of parton energy loss for both light and heavy quarks 
| traversing the hot and dense medium remains illusive. 
^ . It is of great interest to investigate the path length (L) and medium density dependence of the energy loss for partons 
' traversing the hot and dense medium created in nucleus-nucleus collisions. The collisional energy loss per unit length 
depends on the medium density and on the differential cross section weighted by the energy transfer (l2l.ll3l.[l3.[l5lfla|. 
In case of the radiative energy loss, gluons are radiated in medium-induced multiple scattering processes, and the 
energy loss has been calculated to depend, to leading order in 1/E, on L and on the average squared transverse 
momentum transferred to the hard parton per unit length 0, d, EE EH Eli EH Ell ■ Thus both radiative and collisional 
energy losses are closely related to L and the soft parton rapidity density of the medium (x^^-)- If we assume the 
linear dependence on L and -jj- ^jj- , both approximately proportional to the cube root of the number of participants 
(Np ar t), then the energy loss goes with as observed for 7r° in Ref. On the other hand, if is approximated 

by the experimentally measured A± by the overlap area of the collision system S, and L by \/~S, then the energy 
loss goes with ^ . Both N part and S can be estimated with Monte Carlo Glauber calculations [2(j • In this paper, 
we extract the effective fractional energy loss for high px particles (5 < px < 10 GeV/c) using published data from 
STAR and PHENIX experiments [H, 0, HL@> an( i study the centrality dependence. The px region has been selected so 
that we can better focus on the final state medium-induced energy loss since the initial-state effects, such as nuclear 
modification of the parton distribution functions and parton intrinsic px broadening, are expected to have small 
impacts [HI, HI]. 

Nuclear effects on the particle spectra are studied by comparison with a nuclcon-nuclcon (NN) reference via the 
nuclear modification factor [23| : 

d^N AA /dpxd V 

RAAipT) = TAAW/dprdr,' (1) 
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FIG. 1: (color online) STAR charged particle spectra in Au+Au (0 — 5%) and p+p collisions at 200 GeV [j§]. The insert shows 
the high pr region in a log-log scale. The fitting curves are in the form of the power law function, as discussed in the text. 




FIG. 2: (color online) Fractional energy loss Apr/pr (for pr > 5 GeV/c) obtained from Eq. ((5JI versus centrality given by 
N^ t . The Raa values are from Refs. @, H H, 0. The data points for 200 GeV Au+Au (Cu+Cu) collisions are fit with a 
straight solid (dashed) line. The /ndf — 30/19 when we fit all the data points with a single straight line (not shown in the 
figure). 



where Taa = (-^fam)/^^ from a Glauber calculation accounts for the nuclear collision geometry. STAR's mea- 
surements have shown that for pt > 5 GeV/c, Raa is roughly constant and significantly lower than unity, for both 
charged particles 0] and non- photonic electrons Q. The same fact holds for the n° measurements by PHENIX 0, 5]. 

In Fig.[TJwe show a comparison of STAR's charged particle spectrum between Au+Au (0—5%) and p+p collisions 3j, 
and fit the data points for px > 5 GeV/c with the power law function 

-d 2 N/dp T dr, = A(l+p T /p )- n . (2) 



2irpT 

In the insert, fitting functions are straight lines in the log- log plot. The suppression in the nuclear modification factor 
can be viewed as a horizontal shift effect [1] in the spectrum from p+p to Au+Au collisions: 



d 2 N AA ( PT )/d PT d v = d 2 a NN {p' T =p T + S(p T )) dp' T _ d 2 a NN (p' T ) [i | dSjpr) 
Taa dpipdr] dpr dp' T d?] dpx 



(3) 



With Eq. © and Eq. ©, Eq. (TJ) now becomes 



W = g ^j^rt], (4) 
(I+Pt/Po) n PT dp T 



If we assume the fractional px shift in the spectrum is a constant, S{pt)/pt = So, then 



R A Ap T )= { ^ +PT / p :!:\ i+s G r^. ( 5 ) 

(I+PT/PO) " 
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FIG. 3: (color online) Fractional energy loss Apr/ 'pt (for pr > 5 GeV/c) obtained from Eq. <(Sj versus The Raa and 

=P- values are from Refs. 

S 3 0, Si and Refs. 0, HH, respectively. The fit is in the form of a straight solid (dashed) line 
for 200 GeV Au+Au (Cu+Cu) collisions. The X 2 /ndf — 30/19 when we fit all the data points with a single straight line (not 
shown in the figure). 



The parameters po and n are obtained and fixed when we fit the spectrum in p+p collisions as shown in Fig. [TJ so 
the Raa of high px particles is closely related to So- In case of small po and So, 

Raa ~ (1 + SoT n+2 (6) 

is a constant. In reality, the typical value of po is between 0.2 and 0.3 GeV/c, and So is smaller than 0.3. 

The effective fractional energy loss Apr/pr is related to the fractional shift in the measured spectrum, So: 
Apx/pr = So/(l + So). For px > 5 GeV/c in 200 GeV Au+Au collisions, we extract the Apx/pr values by fit- 
ting with Eq. ([5} STAR Raa(pt) results for both charged particles [|[ and non-photonic electrons Q, and PHENIX 
Raa(pt) result for 7r° Figure [2] and [3] show the effective fractional energy loss versus centrality given by Np^t 
and j respectively. The values of ^ in 200 GeV Au+Au collisions arc from Ref. j24j]. The vertical error bars 

represent the statistical errors obtained in the fitting procedure. In both figures, the energy losses of charged parti- 
cles, 7T° and non-photonic electrons, despite different fragmentation functions and decay processes involved for these 
particles, follow the same curve. In parton energy loss scenario, it is generally believed that partons would lose energy 
in the hot and dense medium and the fragmentation process takes place after the partons escape into the vacuum. 
Assuming that the vacuum fragmentation process is the same in p+p and in A+A collisions, different particles would 
reflect the same amount of leading parton energy loss. This scenario is supported empirically by the fact that r\ and 
7T , proton and charged tt have the same Raa values at the high pt region from PHENIX and STAR measurements 
[2a , l2q . The semi-leptonic decay kinematics in the heavy quark decays will not change the fact that the transverse 
momenta of non-photonic electrons are on average proportional to those of heavy quarks. Therefore, our extracted 
Apr Ipt ratios represent the fractional parton energy loss in the medium. 

For Au+Au data with various collision centralities the fractional parton energy loss, Apx/pr, increases approxi- 
mately linearly with Np^ t in Fig. El and with ^ in Fig. [31 These N part and ^ dependences are consistent 

with the scenario of energy loss being determined by path-length-times-density as suggested by GLV (2?| and PQM 
[28j |. However, the accurate Cu+Cu data on 7r° from PHENIX Q covering the low N part region, as shown in Fig. [2] 
and[3j systematically deviate from the suggested dependence, calling for possible alternative path length and medium 
density dependence for parton energy loss. Here, the values of ^ in 200 GeV Cu+Cu collisions are converted from 

PHOBOS measurements of ^ (29|, scaled by a factor of 1.15 from pseudo-rapidity density to rapidity density [30j . 

In Fig. [H we show Apr /pt as a function of the particle rapidity density per transverse area (^^r), which appears 
to be a better quantity to describe the dependence of parton energy loss in the hot and dense medium: the values of 
fractional parton energy loss, Apx/pr, from Au+Au and Cu+Cu collisions seem to depend linearly on ^ =^ within 

statistical errors. The fitting line intercepts at a finite -57^ for zero parton energy loss, possibly indicating that a 
minimal medium density and/or transverse dimension of the colliding system is needed in order to have final state 
energy loss. We note that the observed dependence for parton energy loss implies the absence of or weak 

path length dependence for parton energy loss. The particle rapidity density per transverse area, a quantity likely 
proportional to the initial parton density of the collision, may determine the magnitude of the parton energy loss. 
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FIG. 4: (color online) Fractional energy loss Apt/pt (for pr > 5 GeV/c) obtained from Eq. ((Sj versus centrality given by 
^T^r- The Raa and ^ values are from Refs. 0, 0, H El and Refs. [24l. 129] . respectively. The fit is in the form of a straight 
line for all data points. When we fit the data points in Au+Au or Cu+Cu separately, the fitting results are consistent with the 
result shown in the figure. 

The absence of the path length dependence, if confirmed, would contradict calculations of parton energy loss that 
employ a static average geometry. 

Possible physical explanation for the lack of strong path length dependence lies in the fast expansion of the collision 
system, leading to a rapidly dropping medium density as a function of time. We consider the energy loss as a function 
of the expansion time of the collision system, in a scenario proposed by Ref. [3l| using a two-component dynamical 
model (hydro + jet model) [13] with a fully three-dimensional hydrodynamic model [33[ for the soft physics and pQCD 
jets for the hard physics which are computed via the PYTHIA code |34| . Figure 4 of Ref. [3l[ shows a hydrodynamic 
calculation of the jet quenching rate N j et (r) / N j et (To) for px = 5 GeV/c jets in Au+Au collisions at 200 GeV, where 
Njet(T)/Nj e t (to) is equivalent to the Raa of the px = 5 GeV/c jets at the expansion time r. We apply Eq. ([6]) to 
extract the energy loss information from the Nj et (r) / Nj et (ro) (or Raa) curves, and plot the effective fractional energy 
loss against r in Fig. [5] The dashed line corresponds to the "constant energy loss" assumption, that the energy loss per 
unit length is proportional to the local parton density p(r) in the medium, and the solid line, the more sophisticated 
GLV formula [27] . For both energy loss schemes and both central and mid-central collisions, a large portion of energy 
loss occurs within the first a few fm/c. A high px parton may traverse the full path length along the trajectory inside 
a fast expanding medium, but the energy loss will effectively cease after the first a few fm/c making the total path 
length irrelevant to the magnitude of the energy loss. As a result, the effective path lengths are close to a constant 
for all collision centralities of both Au+Au and Cu+Cu collisions, and the parton energy losses depend only on the 
medium density at the early stage of the collision which may be represented by ^^r- Note that the hydrodynamic 
calculation is only valid after thermalization of the system. The actual effective time interval for parton energy loss 
may be shorter than what the model indicates, because of non-equilibrium dynamics in the early stage. Detailed 
calculations may be needed to address why peripheral Au+Au and Cu+Cu data follow the linear dependence 
where the participant matter system is not much larger than a few fm in radius. 

This physical picture for parton energy loss, if confirmed by more precise data over a broad range of N part , would 
have many implications on theoretical calculations for parton energy loss and on azimuthal angular anisotropy of 
high pt particles. Theoretical calculations using static nuclear geometry with average overlapping participant matter 
will not be able to capture the essence of the limiting time duration for parton energy loss processes. Such static 
calculations, for example, Refs. [H, [36J concluded that in central A+A collisions the surviving high px particles 
tend to come mostly from surface region of the participant matter because partons from the central region would be 
quenched along the full path length. However, our energy loss scenario implies that partons from the central region 
would suffer energy loss in the first a few fm/c and may escape producing high px particles in the final state. We 
argue that the energy loss mechanism is intrinsically a dynamical process and cannot be adequately described by 
static calculations. Theoretical calculations with a three-dimensional hydrodynamic evolution model show that an 
energy loss that grows linear (quadratic) in L in a constant medium, as characteristic of collisional (radiative) energy 
loss, translates into a logarithmic (linear) path length dependence in a medium undergoing longitudinal Bjorkcn 
expansion (T5, l37l. l38l. |39| . In other words, the rapid hydrodynamic expansion weakens the path length dependence of 
the energy loss. This theoretical result is consistent with our empirical conclusion that the rapid dynamical expansion 
may account for the path length independence of the energy loss implied by experimental Raa measurements. 

Elliptic flow (V2) is defined to be the second-order harmonic of the Fourier expansion of particle's A<p distribution 
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FIG. 5: (color online) Fractional energy loss Apx/pr (for 4.5 < pr < 5.5 GeV/c) versus the expansion time of the collision 
system. The Nj e t(r)/Nj e t(To) (or Raa) values are from Ref. [3l[, where the dashed line corresponds to the "constant energy 
loss" assumption, and the solid line, the GLV formula [27[, as discussed in the text, b = 2 (7) fm corresponds to — 5% 
(15 - 25%) most central Au+Au collisions, with N var t « 350 (200). 



[40|, where A<p is the particle's azimuthal angle with respect to that of the reaction plane. The high px (> 5 GeV/c) 
particle vi , has been considered as arising from path length differences for parton energy loss as a function of A</>. Our 
dynamical energy loss scenario would suggest that V2 values for high px particles would not be as large as predictions 
based on static energy loss picture. The Raa of high px 7r° has been observed to depend on the emission angle A</> Q], 
an equivalent measurement to that of V2- The reported angular dependence for px 5 — 8 GeV/c seems to be weaker 
than that for 3-5 GcV/c. We cannot draw a firm conclusion on the dynamic scenario for parton energy loss 
without detailed comparison with theoretical calculations. Significant theoretical and experimental uncertainties arise 
from large soft particle contributions to the measured px range and from the event-plane determination. Both may 
have biases in the angular distributions. The event plane is very often estimated with the azimuthal angle distribution 
of the detected final-state particles, and "in-plane" is defined as the direction where most particles come out. If there 
are some extra azimuthal correlations (not related to the reaction plane orientation) between the particles of interest 
and the ones participating in the event plane determination, then the measured V2 will deviate from the true value, 
which is called non-flow effects including momentum conservation [4l| , long- and short-range two- and many-particle 
correlations due to quantum statistics, resonances, mini and real jet production, etc (40l |. Non-flow effects could easily 
result in more particles in-plane than out-of-plane, even if there are no collective motions like elliptic flow or path 
length dependence of the energy loss in the high px region, especially in the most central and peripheral collisions [12] • 
Often to suppress the non-flow effects, the particles used to estimate the event plane are selected to be one or two units 
of rapidity away from the particles under study. But long range correlations may still exist. Further improvement 
on the event plane determination may use the side- ward deflected spectator neutrons [HI, [H| as carried out in STAR 
measurements [45| . 

Another way to study the path length dependence of the energy loss is through the A<p dependence of the di- 
hadron azimuthal correlation after the subtraction of the elliptic flow modulation. STAR's measurements show that 
the away-side correlation evolves from single- to double-peak with increasing A</> [46| . where the high px trigger 
particles range from 3 to 4 GeV/c. However, one difficulty in this analysis lies in the determination of the elliptic flow 
background, which could be larger than the correlation of interest by two orders of magnitude. Non-flow effects could 
influence both the measured V2 values and the symmetry of the background correlation, and flow fluctuations bias the 
event plane resolution, an important quantity to calculate the flow-induced two-particle azimuthal correlations in this 
analysis 47 1. Further, the Zero Yield at Minimum (ZYAM) approach (48| which has often been used to normalize the 
background introduces additional uncertainties. In addition, the analysis requires large statistics and current results 
are limited to leading particles with 3 < px < 4 GeV/c, which is below the px domain we considered relevant for 
parton energy loss study. Correlation studies triggering on leading particles above 5 GeV/c are needed to test the 
physical scenario we have proposed. 

In summary, we have extracted the effective fractional energy loss Apx/px for high px charged particles, non- 
photonic electrons and ir° in Au+Au collisions and ir° in Cu+Cu collisions at 200 GcV. Empirically the Apx/px is 
found to be a constant for flat nuclear modification factor Raa as a function of px from a given collision centrality. 
The derived fractional parton energy loss Apx jpx seems to depend on particle rapidity density per transverse area for 
Cu+Cu and Au+Au data over the full range of N part , which implies that there is no strong path length dependence 
for parton energy loss along the trajectory of N part geometry. We argue that the absence of strong path length 
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dependence may be due to the rapid expansion of the participant matter in the early stage so that parton energy 
loss only takes place within the first a few fm/c duration. Simultaneous studies of nuclear modification factors Raa 
and elliptic flow v- 2 for high (> 5 GcV/c) particles would shed more insight on this dynamical scenario for parton 
energy loss. 
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